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Abstract. Homogeneous isotropic cosmological models with two torsion functions filled 
with scalar fields and usual gravitating matter are built and investigated in the framework 
of the Poincare gauge theory of gravity. It is shown that by certain restrictions on indefinite 
parameters of gravitational Lagrangian the cosmological equations at asymptotics contain 
5^ ' an effective cosmological constant that can explain observable acceleration of cosmological 

expansion. The behavior of inflationary cosmological solutions at extremely high energy 
densities is analyzed, regular bouncing solutions are obtained. The role of the space-time 
04 ' torsion provoking the acceleration of cosmological expansion is discussed. 
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, 1. Introduction 

X 

d • The discovery of the acceleration of cosmological expansion at present epoch is the most 
principal achievement of observational cosmology at last time [1]. By using Friedmann 
cosmological equations of General Relativity theory (GR) in order to explain accelerating 
cosmological expansion, the notion of dark energy (or quintessence) was introduced in 
cosmology. According to obtained estimations, approximately 70% of energy in our Universe 
is related with some hypothetical form of gravitating matter with negative pressure — "dark 
energy" — of unknown nature. Previously a number of investigations devoted to dark 
energy problem were carried out (see review [2]). According to widely known opinion, the 
dark energy is associated with cosmological term. If the cosmological term is related to the 
vacuum energy density, it is necessary to explain, why it has the value close to critical energy 
density at present epoch (see for example [3]). 
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The present paper is devoted to investigation of the "dark energy" problem in the 
framework of the Poincare gauge theory of gravity (PGTG), which is a natural generalization 
of Einsteinian GR by applying the gauge approach to the theory of gravitational interaction 
(see for example [11]). In fact the generalization of GR leading to PGTG is necessary, 
if one supposes that the Lorentz gauge field corresponding to fundamental group in 
physics - the Lorentz group - exists in the nature (see [4]). According to PGTG the 
physical space-time possesses the structure of Riemann-Cartan continuum with curvature 
and torsion. The investigation of isotropic cosmology built in the framework of PGTG (see 
[4-6] and references herein) shows that the gravitational interaction in PGTG, unlike GR 
and Newton's gravitation theory, can have the repulsion character in the case of gravitating 
systems with positive energy densities and pressures satisfying energy dominance condition. 
So, the gravitational repulsion effect takes place at extreme conditions (extremely high 
energy densities and pressures) preventing the appearance of cosmological singularity in 
homogeneous isotropic models (HIM) [5]. According to generalized cosmological Friedmann 
equations (GGFE) for HIM deduced in the framework of PGTG, all cosmological solutions 
including inflationary solutions are regular in metrics, Hubble parameter, its time derivative 
and have bouncing character. Properties of discussed HIM in PGTG coincide practically 
with that of GR at sufficiently small energy densities, which are much less in comparison 
with limiting (maximum) energy density for such models. By including cosmological term of 
corresponding value to GCFE, we can obtain regular cosmological solutions with observable 
accelerating expansion stage. However, like GR, the problem of dark energy in such theory 
is not solved. 

From geometrical point of view, the structure of HIM in PGTG can be more complicated 
in comparison with models describing by GGFE. In fact in the case of homogeneous isotropic 
models the torsion tensor S^^i, — —S^^^ can have the following non-vanishing components 

[7, 8]: S^io = S^2o = S^3o = Si{t), 5*123 = 5*231 = 5*312 = 5*2 (t) ==sin6', where 5*1 and 

Vl — A'r^ 

5*2 are two torsion functions of time, spatial spherical coordinates arc used. The functions 5*1 
and 5*2 have different properties with respect to transformations of spatial inversions, namely, 
the function 5*2 (t) has pseudoscalar character. The GGFE follow from gravitational equations 
of PGTG for HIM together with 5*2 = 0. Obtained physical consequences of GGFE have 
principal character. However, it is necessary to note that gravitational equations of PGTG 
for HIM have also other solution with non- vanishing function 5*2. 

The HIM with two torsion functions filled with scalar fields and usual gravitating matter 
are studied below in the frame of PGTG in connection with the dark energy problem. 
Following to [9] , in Section 2 cosmological equations for such HIM are introduced. In Section 
3 the solutions asymptotics of cosmological equations is analyzed. In Section 4 the bouncing 
character of inflationary cosmological solutions is examined. 
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2. Cosmological equations for HIM with two torsion functions in PGTG 

At first, let us mention some general relations of the PGTG. Gravitational field is described 
in the frame of PGTG by means of the orthonormalized tetrad h\ and anholonomic Lorentz 
connection A'^^^ (tetrad and holonomic indices are denoted by latin and greek respectively); 
corresponding field strengths are torsion S'^ni, and curvature F'^'^^^ tensors defined as 

We will consider the PGTG based on the following general form of gravitational Lagrangian 

CG = h[foF + F^f^^^ (/i + h Fa,0, + h F^.o.p) + F^^^ {f, F^, + f, F,,) + f, F^ 

+3'^^" (a, S^,, + a2 S,,^) + S'\,^Ss^^] , (1) 

where h = det {h\), F^^ = F'^^au, F = F^^, fi {i = 1,2, ... ,6), ak {k = 1, 2, 3) are indefinite 
parameters, /o = (IGttG)"^, G is Newton's gravitational constant. Gravitational equations 
of PGTG obtained from the action integral I = J {Cg + Cm) d'^x, where Cm is the Lagrangian 
of matter, contain the system of 16+24 equations corresponding to gravitational variables 
h\ and 

Any homogeneous isotropic gravitating system in PGTG is characterized in general case 
by three functions of time: the scale factor of Robertson- Walker metrics R and two torsion 
functions and 5*2. Below the spherical coordinate system is used and the tetrad is taken in 
diagonal form. Then the curvature tensor has the following non- vanishing tetrad components 
denoted by means of the sign ' : 

p6i _ p02 _ pQS _ A 77.12 _ jpld, _ rp23 — A 

^ 6i ~ -'^ 02 ~ -'^ 03 — ^1' ^ 12 — ■'^ 13 ~ ^ 23 — ^2, 

pOi,, _ p62 z^OS^ ^ A p32 c^i3 p2i_ = A 

^ 23 ~ ^ 31 ~ ^ 12 — ^3) 01 — ■'^ 02 ~ ^ 03 — ^4 

with 

Ai^H + H^-2HSi-2Si, 

^2 = + (i/ - 2S,f - Si ^2) 

As^2{H-2Si)S2, 

A4 — S2 + HS2, 

where H — R/R is the Hubble parameter and a dot denotes the differentiation with respect 
to time. 

The Bianchi identities in this case are reduced to two following relations: 
A2 + 2H {A2 - Ai) + ASiAi + 2S2A4 = 0, 
ia + 2// (^3 - A4) + 4,^1^4 - 2,^2^ = 0. 
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The system of gravitational equations for HIM corresponding to gravitational Lagrangian 
([1]) has the following form 

a{H- Si) Si - 2bSl - 2/0A2 + 4/ [Al - Al) + 2q2 {Al - Al) = 



a{Si + 2HSi - Sf 



2hSl 



2/0 {2A, + A2) - 4/ {Al - - 2g2 [Al - A^) = p, 



/ 



(I2 



ii + 2H (Al - A2) + ASiA^ 



^25*2^3 - gi 5*2^4 + /( 



5i = 0, 



i4 + {A, - A3) + 451^3 - 4/ ^2^2 - 2gi52Ai - (/o - 6) ^2 = 0, 



(4) 
(5) 
(6) 
(7) 



where 



a = 2a\ + 02 + 803, 

h 



02 ~ 0,1, 



/=/l + y + /3 + /4 + /5 + 3/6, 

gi = /2 - 2/3 + /4 + /5 + 6/6, g2 = 2/1 - /2, 

p is the energy density, p is the pressure and the average of spin distribution of gravitating 
matter is supposed to be equal to zero. Equations (!lj)-(l5l) lead to generalization of Friedmann 
cosmological equations of GR, which does not contain high derivatives for the scale factor i?, 
if a = (see below). Moreover, equations ([6])-([7]) take more symmetric form, if 2/ = gi + g2- 
Then by using the Bianchi identities (JSj), the system of gravitational equations for HIM take 
the following form: 

1 



- 26 Si - 2/0A2 + 4/ [A\ - A^) + 2g2 [A\ - A^) : 

- 26 ^2^ - 2/0 (2Ai + A2) - 4/ [A\ - A^) - 2g2 [Al - A^) = p. 



/ 



ii + i2 + 45i (Al + A2) + g2 ^2 (^3 + Ai) + /o5i = 0, 



(9) 
(10) 



g2 [(A3 + A4) +4Si(A3 + A4)] - 4/ S2(Ai + A2)-(/o- 6)^2 = 0. (11) 

The system of equations (l8l)-(|TTl) together with definition of curvature functions ((21) is the 
base of our investigation of HIM below. Note also that the conservation law for spinless 
matter has the usual form: 

p + 3i/(p + p) = 0. (12) 

In order to investigate inflationary cosmological models we will consider below HIM filled 
with non-interacting scalar field minimally coupled with gravitation and gravitating matter 
with equation of state in the form -pm = Pm{pm) (values of gravitating matter are denoted 
by means of index "m"). Then the energy density p and the pressure p take the form 



(13) 
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where V = V{(f)) is a scalar field potential. By using the scalar field equation in homogeneous 
isotropic space 

dV 

<P + ^H<P=-— (14) 

we obtain from f|T2|) - f|T3|) the conservation law for gravitating matter 

+ 3/7 (pn, + = 0. (15) 
From (IH])-© follows that 

A^ + A, = -^{p-?,p)-^Sl. (16) 
J-^/o Jo 

By using ( ITGil and the formula following from definition of curvature functions and ^4 



/I3 /I4 



i?2 



we find from gravitational equations 



5*2 + HS2 



the following expressions for Ai and A2: 



Ai 



1 



12/0^ 



p + 3p--{p~3p-12bSi) 



ae 



If 

3aefo 



HS2 + S2] +4 



i?2 



- s ' si 



(17) 



A, 



6foZ 



P - Qbsl + ^ {p - 3p - mslY 



3aefo 



HS2 + S2] +4 



k 



- S' S 



where Z = 1 + a {p - 3p - 12 {b + efo) S|) = 1 + a (AV - (j)^ + p^ - 3p^ - 12 (6 + efo) Sj j , 
a = //(3/q), e = ^2// (hence, q2 = 3a ef^). By using (fT^ - (ITBI) and the following relation 
obtained from definition of ^3 and A4 

As + A^ = S2 + 3HS2 - 451^2, (19) 
we find for the torsion function Si from (fTOj) the following expression: 



Si = - 



3a 



dV ■ 



where 



Y = {p^+p^){3 



+ H{Y + 

dPm 
dpm 



-4(26-5/0)^2^2 



(20) 



1 + 12£/o^2- 



Then by using formulas flTBl) and flTI?]) we find from (ITTl) the following second order differential 
equation for the torsion function 5*2 : 



^2 + 3^^^2 + 3HS2 - 4 5i - 3HSi + 4^1^ S; 



^^{p-3p-mSl)S2-^-^^S2 
-J Jo / 



0. 



(21) 
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The obtained expressions (fT7l) - (ITg|) for curvature functions A2 and Ai together with their 
definition ([2]) give the generahzation of cosmological Friedmann equations for HIM: 



i?2 



+ {H~ 2Si) 



1 



QfoZ 
3aefo 



p + 6 ifoZ - b) Si + ^ (p - 3p - UbSiy 



HSo + So] +4 



H + H^- 2HSi - 2Si 



l2foZ 



^ q2 



3p - ^ {p - 3p ~ msiy 



Si 



(22) 



{p - 3p - 12bSi) Si + ^ 



2 

Saefo 



HSo + So] +4 



k 



-Si 



S' 



(23) 



These equations contain the torsion function 5*1 determined by fl20l) and the torsion function 
S2, satisfying the equation fl^ . Obtained equations contain three indefinite parameters: 
indefinite parameter a determining the scale of extremely high energy densities [4] , parameter 
b with dimension of parameter /q and the parameter e without dimension. We have to 
analyze all these equations in order to investigate HIM with pseudoscalar torsion function 
in the frame of PGTG. 

3. Asymptotics of cosmological solutions for HIM with pseudoscalar torsion 
function 



The structure of obtained equations fl20l) - fl23l) describing HIM with two torsion functions is 
essentially more complicated in comparison with the case of HIM with vanishing function 
5*2 . Note that if 5*2 = the equation ( 1211) vanishes and the cosmological equations ( I22i) -( l23l) 
are transformed to GCFE containing the only indefinite parameter a [4,5]. 

Now we will analyze the following question: by what restrictions on indefinite parameters 
the cosmological solutions for HIM with pseudoscalar torsion function have the asymptotics 
in agreement with actual observations. By taking into account that various parameters of 
HIM have to be small at asymptotics, when values of energy density are sufficiently small, 
we see from (1211) . that if \e\ <^ 1, the pseudoscalar torsion function has at asymptotics the 
following value: 

fo{fo-b) P-Sp 



S' 



(24) 



4/6 126 

Then we have at asymptotics: Z (6//0), Si ^ and the cosmological equations (l22l) - 
f l23|) at asymptotics take the form of cosmological Friedmann equations with cosmological 
constant: 



1 

66 



P 



3 (/o - by 

4/ 



(25) 
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H + H' 



12b 



p + 3p 



3 (/o - by 
2/ 



(26) 



From equations fl25|) - fl26l) we see, that parameter b has to be very close to /o, but smaller 
than /q. The value of b leading to observable acceleration of cosmo logical expansion depends 
on the scale of extremely high energy density defined by a~^. If we take into account that 
the value of |(/o — &)^// = \a^^{l — b/ fo)"^ is equal approximately to 0.7pcr (the critical 
energy density is pcr = ^IoHq, where Hq is the value of the Hubble parameter at present 
epoch), then we obtain that b = [1 — (2.8pcra)^^^]/o- If we suppose that the scale of extremely 
high energy densities is larger than the energy density for quark-gluon matter, but less than 
the Planckian energy density, then we obtain the corresponding estimation for b, which 
is very close to /q. Obtained restrictions on indefinite parameters will be used below for 
investigation of inflationary cosmological solutions. 



4. Regular inflationary cosmological solutions with two torsion functions 



To obtain cosmological solution by integrating cosmological equations we have to use the 
equation of state of gravitating matter, which is different at different stages of cosmological 
evolution. So, at asymptotics one uses usually equation of state for dust matter (p = 0). At 
the same time, in order to obtain cosmological solution for inflationary HIM, we will use at 
the beginning of cosmological expansion the expressions (fT3|l for energy density and pressure 
by including scalar field as one component of gravitating matter. Like GR, the inflationary 
stage appears, if the value of scalar fields at the beginning of cosmological expansion is 
sufficiently large (</>>! Mp, Mp is Planckian mass) [10]. 

In order to investigate inflationary cosmological solutions at extremely high energy 
densities, by using (fT3!) - (fT5l) and (l20l) we transform cosmological equations (!22|) - (!23|) and 
equation (|2T]) for 5'2-function to the following form 

1 2 



Y + 



+6aH 
1 



Z + -a[Y + 2(t? 



dV ■ 



^fo) S2S2 



^fo) S2S2 



dV ■ 

Pm + ^<P + V- QbSl + 



^2^ + 4 



k 



Si 



Si 



+ AV -<f- 12bS 



(27) 
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H 



9a 
_3a^ 



' dpi 



2Z 



{Pm+Pmf + 3 



9a2 
2Z2 



-11 



Y + 202) (y + 20' - 12efoSfj 



dpn 
dpn 



(Pm+Pm) +40^ 



4^0 + 2 (26 - 7efo) S2S2 



+ 



3a 



(^|^0-2(2fe-£/o)52^2 



X (r + 202 - 12efoS2S2) + (y + 20^) x (^|^0 - 4 (6 + s/o) 52-^2^ | 



3q; d^V 



902 



dV 



-2(26-5/o) (5| + '52^2)} 

Pm + 3pni - 2 - (/ 



12/0^ 



1 



- -O Pm - ^Pra + 4^ - 0' - 126^2 



+3 



Z 

aefo 



Pm - 3p„, + 4^-02- 126^1 5| 



(28) 



l-^{2b-eU)Si 



+ 3HS2 



l.|(K.2,^)] 



-9^H^S2 



F + 602 + 3^(p^+p^) 



+ 3 



, dpn 
dpn 



-11 + 



a 



CtPrr, / ^ 



+ 3HS2 <! 1 - 4^ f 4^0 - 3 (26 + efo) S2S2 



x(y + 202 _ i2sfoS, 

(^^0 - 2 (26 - £/o) ^2^) (y + 202 - 12£/o52-^2) 

+ (y + 202) (^1^0 - 4 (6 + £/o) ^2^2) } 



-9|,.2 



-e^S2 

z 



H{Y + 2cj>')+2 (^1^0 - 2 (26 - efo) S2S2 
^' + 2 {2b -efo) Si 



dV 
90 



902 
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6a fdV ■ 
+Y (l^*^ - 2 (26 - efo) S,S, 



dV ■ 



3/o 



Pm - + 41^ - 



126^2^ ) + 



-4(6 + 5/0)^252 

fo-h 



f 



52 = 0. 



Equation fl27j) can be written as 

AH"^ + 2BH + C = 0, 

where 



(29) 



(30) 



A 



B = 3a 



C = 9a' 



+ 3q;£:/o5|Z, 
dV ■ 



Z + -a[Y + 2(t? 



^0 - 2 (26 - efo) 
1 , 



2 (26 - 2£/o) 52^2 



+ efQS2S2Z 



+ 3a£:/o 



5,^ + 41-^-5^5. 



1 

'6^ 



Pm + ^0' + - 6652 + _^ (^^ _ + 4r - 



1265^ 



+ 



^ r.2 



From (1301) we obtain two i/±-solutions for the Hubble parameter 



A 



(31) 



where D = B"^ - AAC. 

At asymptotics //_-solutions and iJ+-solutions describe the stages of cosmological 
compression and expansion respectively [4]. The transition from if_-solution to if+-solution 
takes place when D = 

Now we will analyze extremum surfaces in space of independent variables 0, 0, 52, 52, 
Pm, in the points of which the Hubble parameter vanishes H = 0. Extremum surfaces depend 
on indefinite parameters a, e, b and in the case of open and closed models also on the scale 
factor R. By denoting values of variables on extremum surfaces by means of index 0, we 
obtain from (!30|) the following equation for such surfaces 



1 

6^ 



1 • 1 

PmO + 2*^0 + ^0 - 66520 + -a ( Pmo - 3pmo + 4VJ; 



- 0^ - 1265^0 



-9a^ 
-Saefo 



2(26-£/o)52o52o 



k 



'^2 + 4 I ^ - 54 I 5; 



Rl 



20 ^^20 



_ _ Q2 
JTLq 



20 ^0 



(32) 
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where 



Zn = 1 + a 



PmO - 3pmO + 4^0 



12{b + efo)S, 



The derivative of the Hubble parameter on extremum surfaces obtained from (l28l) - (!29l) 
is determined as 



PmO 



3PmO + 2 ' ^ 



1 + a (p 



-'mO 



3pmO + 4Vo 



1 (PmO + Pmo) + 4K) + 



2/o 



2 (PmO - Pmo) + Vq 



-46S'2o + |« ( PmO - 3pmO + 4Vo 



1265^0 



+3aZo 



I PmO 



3pmO + 4K) 



+ (46 - 3£/o) - 4e/o 
2 



i?0 



r.2 



r,2 
•^20 



r.2 
'-'20 



^ r.2 



20 



■-0 



72 
^0 



-108a'£/o52o52o 



— 0o-2(2&-£/o)52o52O 



+6a (26 - efo) S^q \ 72a 



2 (26 - efo) ^20^; 



20 



9V\ • • ■ 

— 00 - (46 + efo) S20S20 
9<l)Jo 



+2 i2b~efo)S',, 



+ 6a Zo 



\d<t>J, 



2 



3/0 



PmO - 3pmO + ^Vo - 00 - 1265*20 + 



/o 



/ 



f33) 



In the case of HIM without pseudoscalar torsion function 5*2 the equation fl32|) and the 
formula ([33]) simplify and take the form obtained in ref. [4]. As it was noted in [4,6], in 
this case the most part of extremum surfaces play the role of bounce surfaces (i^o > 0) 
for scalar field potentials applying in theory of chaotic inflation. The different situation 
is in considering case of extremum surfaces (l32l) . By given values of parameters a and ^ 
the bounce {Ho > 0) takes place only in limited domain of extremum surfaces fl32l) with 
negligibly small values of 820- In the case S20 = the equation of extremum surface ( l32i) 
and the expression (!33|) of derivative Hq are simplifled and take the following form 

If 1 • 1 / • N 2" 

^ PmO + -<^0 + + -a PraO " Sp^o + 4^) " 00 

6/0 L 2 4 V 



I According to the conclusion obtained in Section 3 the parameter b is very close to /o- As result we put 
below for numerical calculations & = /o • 
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/ 



72 



— 3Q;e/oS'2o-^o, 



11 

(34) 



1 

Wo 



(PmO - P: 



+ Vo + -a 



+3a 



dv_y 



1 + a 



PmO — 3pmO + 



PmO - 3pmO + 4\/o - 

4 + (46 - 3£/o) S|o 




(35) 



1 (PmO + Pmo) + 4Fo + 



We see from ( l35l) that the presence of 5*20 in this expression does not prevent from the bounce 
reahzation. Moreover, if we put = and k = 0, from follows that e > 0. 

As an example of inflationary cosmological solutions we will consider below fiat HIM 
filled with ultrarelativistic matter = and scalar field with quadratic potential 

V = |m^0^. For numerical calculations we will use m = 10~^Mp and = 1.2 x 10^^'^ Mp. 
To perform the numerical integration of equations f |T^ -(|T5l). ( I28l) -(l29l) it is convenient 
to transform all variables and parameters entering these equations to dimensionless units 
marked by tilde 



t - 
P - 

H 



t = t/ ^/f^, 

' p = ap, 



R- 

p 

m - 

Sl,2 



R 



p = ap, 



(36) 



The explicit form of equations f|T^ - f|T5l) . fl28|) - fl29|) after this transformation is similar to 
original form except the fact that parameters a and /o are cancelled in obtained equations. 
Particular numerical solution was found under the following value of indefinite parameter 
e = 10~'^. Initial conditions for H, (p, S2, S'2, Pm were taken at a bounce as follows 



Ho = 0, 



25, ^20 = 0, S' = 0.001, 



PmO 



0.4, 



where the prime denotes the differentiation with respect to t. Initial condition for 0q was 



taken to satisfy fl51|) . Obtained solution is given in figure 1 - figure S] and includes four 
stages: the compression stage (figured]), the transition stage from compression to expansion 
(figure [2]), the inflationary stage (figure [3]) and the postinfiationary stage (figure H]). The 
distinguishing features of obtained solution are its completely regular character. Note, that 
during inflationary stage number of e-folds for the scale factor is equal approximately to 76. 

Similar to GR, the transition to radiation dominated stage can be realized by 
transformation of oscillating scalar fields (see figure 4) into particles [11]. Details of such 
transition in considered theory require further investigation. In particular, the presence of the 
oscillations of the Hubble parameter H (figure 4) can lead to some distinguishing features 




Figure 1. Compression stage. 



of the inflationary scenario in the considered theory in comparison with the inflationary 
scenario in GR. 

After transition to matter dominated stage the further evolution of the Universe in this 
theory is the same as in the frame of standard cosmological scenario. The transition to 
the accelerating expansion takes place, when the value of effective cosmological constant is 
greater than the matter energy density. 



5. Conclusion 



As it was shown, in the framework of PGTG the gravitational interaction in the case of usual 
gravitating systems can have the repulsion character not only at extreme conditions [4,5], but 
also at sufficiently small energy densities. The pseudoscalar torsion function in HIM provokes 
the appearance of effective cosmological constant at asymptotics of cosmological solutions 
that can lead to observable accelerating cosmological expansion. Quantitative agreement 
of the obtained result with observations depends on corresponding restrictions on indefinite 
parameters a, b and e from Section 3. Numerical solution for inflationary cosmological model 
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Figure 2. Transition stage. 

presented at Figures 1-4 conserves its qualitative behaviour by relatively small variations of 
indefinite parameters and initial conditions. 

The effect of acceleration of cosmological expansion in PGTG has the geometrical nature 
and is connected with geometrical structure of physical space-time. Hence, from the point 
of view of considered theory hypothetical form of gravitating matter — dark energy — is 
fiction. 
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Figure 3. Inflationary stage. 
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Figure 4. Postinflationary stage. 



